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SUMMARY 


Integrating aatrices for* the basis of an efficient numerical proce- 
dure for solving differential equations associated with rotating bea« 
configurations. In vibration probleas, by expressing the equations of 
notion of the bean in natrix notation^ utilizing the integrating matrix 
as an operator, and applying the boundary conditions, the spatial depend- 
ence is removed from the governing partial differential equations and 
the resulting ordinary differential equations can be cast into standard 
eigenvalue fora. Previous derivations of integrating matrices have been 
restricted to one spatial dimension. This report derives integrating 
matrices based on two-dimensional rectangular grids with arbitrary grid 
spacings allowed in one direction. The derivation of higher dimensional 
integrating matrices is the initial step in the generalization of the 
integrating matrix methodology to vibration and stability problems 
involving plates and shells. 


INTRODUCTION 

The equations of motion governing the vibration and stability of 
rotating structures have no closed-form solutions, and approximate 
methods of analysis such as as>'mptotic techniques, finite-element 
methods, or direct numerical integration must be employed to effect a 
solution. For the special case of beamlike (i.e. one-dimensional) struc- 
tures. a numerical procedure based on the use of "one-dimensional" 
integrating m.itriccs (.refs. 1, 2) provides a particularly efficient 
procedure for solution (refs. 3-5). An integrating matrix provides a 
means for inte^r.»ting a function whose values arc known on a grid composed 
of increments in the independent variable. Such matrices can be derived 
by approximating the integrand by an appropriate polynomial (refs. 2, 5). 
In bc.ara vibration problems, for example, expressing the equations of 
motion in matri' notation, utilizing the integrating matrix as an operator 
and applying the boundary conditions remove dependence on the single 
spatial variable, and the resulting ordinary differential equations can 



be put into standard eigenvalue form. Solutions can then be determined by 
standard methods. 

The integrating matrix technique for numerically integrating a func- 
tion defined by its values at a set of grid points separates data defining 
the function at the grid points from data defining the structure of the 
grid. For one-dimensional integrating matrices, i.e. matrices for inte- 
gration with respect to a single space variable, data defining the grid 
structure uniquely deteraine the elements of the matrix. Thus, as long 
as the same grid points and same degree approximating polynomials are used, 
the same integrating matrix may be used to integrate any number of dif- 
ferent functions. This property allows the integrating matrix to be used 
as an operator. This is in contrast to numerical integration schemes 
based on interpolation or divided differences where the analogous quad- 
rature formulas involve a combination of values of the function at the 
grid points as well as data defining the structure of the grid. 

Although the derivation and use of one-dimensional integrating 
matrices has been treated in the literature (refs. 1-5), the derivation of 
a "two-dimensional" integrating matrix for integrating a function of two 
variables whose values are known on a two-dimensional grid has apparently 
not been addressed. The purpose of this report is to derive integrating 
matrices based on two-dimensional rectangular grids with uniform spacing 
in one of the two orthogonal directions. As in the one -dimensional 
integrating matrix, the elements of the two-dimensional integrating matrix 
depend only on the data defining the grid structure and the approximating 
polynomials. 

PRELIMINARY CONSIDERATIONS 

Derivation jf a two-dimensional integrating matrix is related to 
numerical approximation of the imiltiple integral 

f(x,y) dx dy (1) 

R 
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where R 


is the rectangular region 


R 



a < X < b. 


c 



and f(x,y) is a continuous function on R. Values of f(x,y) are 
assumed known at the NM (= N M) grid points 


G 


i * 1, .... N; j = 1, .... M 
(x^, y^); a » XI < X: < ... < Xj^ = b 
c = yi < y: < < yj^^ = d 


( 2 ) 


which lie at the intersection of lines dividing the region R into 
(N - 1) (M - 1) rectangles. A portion of the grid is shown in figure 
1. Grid spacing nay be nonuniform in either the x-direction or y-direction, 
but not both. In the present derivation, as shown in figure 1, non- 
uniformity is allowed in the x-direction. 

From the calculus, the double integral [eq. (1)] may be written as 
an iterated integral 

A = I f(x,y) dxjdy 
Hence, if 

F(y) = J**’ f(x,y) dx (4) 

a 


then 


A 



F(y) dy 


(5) 


The variable y i.s held fixed while performing the x-integration in equation 
(4). Equations (4) and (5) thus effectively split the multiple integral 
into two single integrals of functions of one variable. This decomposition 


3 



Figure 1. Portion of the grid G which divides the region R into (N - 
'< (M - 1) rectangles. In this sample grid, the spacing is non 
uniform in the x-direction and uniform in the y-direction with 
y. , = y. + Ay where Ay = (d - c)/(M - 1). 




suggests that two one-dimensional integrating matrices based on the x and 
y grids 

G^: a = xi < X 2 < . . . < = b (6a) 

c = Xi < X2 < . . . < y^ = d (6b) 

might be used to produce a two-dimensional integrating matrix based on the 
rectangular grid G. Before exploring this possibility, it is necessary to 
develop some notation. 

Notation for Data and Integrals 
To describe data defining f(x,y) at the grid points, let 

f.. = f(x., y.) i = 1 N; j = 1 M (7) 

Integrals with respect to x from a = xj to Xj^ for fixed y will now be 
denoted by 

X, 

Fj.(y) =f if(x.y) dx k = 1 N (8) 

xi 

These integrals play the role of F(y) in equation (4). Similarly, the 
role of A in equation (5) will be played by NM integrals of the form 

A - 

yi ‘ly j = 1 , M; k = 1 . N (9) 


where yi = c. Combining equations (8) and (9) gives 



y . X 

^ f ^ f(x,y) dx dy 

yi 


( 10 ) 


The one-dimensional N by N integrating matrix based on the x grid 

will be denoted by [I], while the one-dimensional M by M integrating 
matrix based on the y grid will be denoted by [J]. 
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To approximate integrals of f(x,y) with respect to x for fixed 
y, the N by N matrix [I] premultiplies the N by 1 column vector 
{f^} which, for the fixed y value X * say, contains data defining 
f(x, y^^) on G^. The data required by [I] is shown in figure 2. In 
particular, if the column vector is defined by 

tfj} = (fit f2t ••• fNfe)^ (11) 


and 


(F^} =((0 F2(y^) F3(y^) ... F^(y^)f 


( 12 ) 


then 


[I] ffj = (F.) (13) 

A- X. 

To approximate integrals of Fj^(y) with respect to y for fixed k, 

the M by M matrix [J] premultiplies the M by 1 column vector {H,^} 

which contains data defining F. (y) on the grid G . The data required 

K y 

by [J] is shown in figure 3. For fixed k = 1 N, defining the 

column vector {Hj^} by 

'"k* ■ (fk')"'’ fk'>'2> 

gives 

[J] = (0 Ajk A3k . . . (15) 

A comparison of figures 2 and 3, or equivalently comparing equation 
(14) with the right-hand side of equation (13), shows that the output 
from approximating an x integration is not, by itself, suitable input 
for approximating a y integration. To obtain the data required to 
form any of the vectors x integrations must be performed at all 

M of the y levels, and the resulting values must be rearranged from 
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Figure 2. A grid G with N * 6 and M = 4. Grid points marked with # 

denote points at which data is required for computation of the 
x-integrals FjcCXj) (k = 1, . . . , 6) using the integrating matrix 
[I] based on the grid G . 




Figure 3. A grid G with N = 6 and M = 4. Grid points marked with • 
denote points at which data is required for integrating the 
function F)^(y) with respect to y using the integrating matrix 
[J] based on the grid G . 




grouping by grid row to grouping by grid column, ti.at is, data must be 
rearranged from variable x-fixed y to fixed x-variable y. This places 
a restriction on the format for presenting data defining f(x,y) at the 
grid points. In particular, to derive a two-dimensional integrating 
matrix on the rectangular grid G, the data for f(x,y) cannot be 
given as an M by N matrix. This result was unexpected as the most "natural" 
way of representing data on an M by N rectangular grid is an M by N 
matrix. While sequential application of [I] and [J] to an M by N 
matrix of data does produce an M by N matrix [A] whose elements ap- 
proximate A., in equation (10), the resulting matrix [A] is not an 
J 

integrating matrix. As this deviation of [A] provides motivation for 
choosing the proper format of data defining f(x,y) at the grid points, 
a brief account of its deviation will be given below. 

A Natural Derivation Which Does Not Yield 
an Integrating Matrix 
T 

Let [D] be the N by M matrix of data defining f(x,y) on G, i.e. 

T T 

the ijth entry of [D] is f^^ . The £th column in the matrix [D] 

is thus the column vector {f^} (? = 1 ,..., M) in equation (11), i.e. 

values of f(x,y) on for fixed y = y^^. Hence, by equation (13), 

the 5,th column in the N by M matrix [E] defined by 

[E] = [I] [D]”^ (16) 

is the column vector {F^} in equation (12) . 

The matrix [E] contains data for the functions Ej^(y), k = 1 ,..., N 

(i.e. x-integrals) on the grid G^. As previously discussed, this data 

must be rearranged into a form appropriate for y- integration using the 

integrating matrix [J]. In particular, the kth column of the M by N 
T 

matrix [E] is the column vector fHj^} in equation (14) (k = 1 ,..., N) . 
Consequently, if [Aj is as previously defined, then 

[A] = [J] [E]^ (17) 

Combining equation (17) with equation (16) now gives 
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(18) 


[A] . [J] [D] [II'^ 

While equation (18) uses the two one-dimensional integrating matrices 
[I] and [J] to obtain [A], it is clear that the right-hand side 
of equation (18) cannot be recast to give the result desired, that is, 

[A] » [P] [D] (19) 

which would allow [P] to be interpreted as a two-dimensional integrating 
matrix. 

The above argument shows that data defining f(x,y) on G cannot 

T 

be expressed as an N by M matrix [D] if the desired end product of 
the procedure approximating the multiple integral [eq. (1)] is to be a 
two-dimensional integrating matrix. In the next section, it will be 
shown that a column vector format for the data defining f(x,y) on G 
will lead to a two-dimensional integrating matrix. 

THEORY FOR TWO-DIMENSIONAL INTEGRATING MATRICES 


Results from previous sections show that representing data defining 
f(x,y) on G as an N by M matrix does not lead to an approximation to 
the double integral [eq. (1)] which can be interpreted as an integrating 
matrix times the matrix of function data. In this section an alternate 
approach is considered where data defining f(x,y) on G is represented 
in column vector form. 


Let (f) be the NM by 1 column vector of data obtained by "stacking" 
the M, N by 1 column vectors which by equation (11) give data 

for fixed values of y = y^, (e = 1 , . . . , M) , that is, 

(f) * ({fj}, {fj}, ..., {fj^^})'^ 

“ (fii» ^ri» •••> ^N 1» ^12* •••> ^Nr. ’ 


IM’ 


’ ^NM^ 


( 20 ) 


10 



Similarly, let {H} be the NM by 1 column vector obtained ' y 
"stacking" the N, M by 1 column vectors {H } in equation (14), that 
is, 

(H) ' ({Hi), {H2> 

~ ^^i(yi)» •••» ♦ ^2(yi)» •••» F 2 (yj^j) » 

pN^yi). .... Fn^V)^ 

Note that each vector (f^} in the stacked vector {f} gives values of 
f(x,y) at the fixed values y = y^ and, hence, is exactly the input 
necessary to approximate the N x- integrals Fj^(yj^) = i N) . 

Similarly, each stacked vector {H } in {H} for fixed m contains 

IB 

data appropriate for approximating integrals of F^(y) with respect to 
y. As seen previously, going from {f} to {H} requires not only 
integrations with respect to x for each fi. 'd y level, but also re- 
arrangement of the resulting data Fj^Cy^) from a row format to a column 
format. In the present devivation, the integration with respect to x 
for fixed y and the data rearrangement will be accomplished simultaneously 
by multiplying {f} by an appropriate NM by NM matrix [T] so that 

{•0 = [T] {f} (22) 

A key observation here is that Fj^ly^j^) may be obtained from the N by N 
integrating matrix [I] based on the grid and the column vector 

by taking the usual dot product of the kth row vector of 

[I] with {f^}: that is, 

^ 23 ) 

Xl 

Hence, each vector {H ) in {H} has the form 

m 

<V ■ CV • 'fi'- "J • I’J 

• (24) 
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(m * 1 N), i,e., only the mth row of [I] enters here. Equation 

(24) and the stacked nature of the vectors (f) and {H} now show 
that the required NM by NM matrix [T] in equation (22) may itself be 
written as a ’’stack" of N, M by NM matrices [T^] (s = 1 N) 

where the nonzero portion of each row of [T^] is simply the sth row 
vector II^J of the N by N matrix [I]. Symbolically, we have that, 
if (OJ denotes the 1 by N zero-row vector, then 


[fj 


a 


UJ loj lOJ 

loj II I loj ... 10 ] 

s 5 V ^ 

• : 

10] 

^ N. 

10] -10] "UJ 


(25) 


and 


’[Til 

[T ) 


[T] = 


[Tn] 


(26) 


It is worth noting at this stage that as (Ij] = [0], [Tj] is simply 
the M by NM zero matrix. Also, in the stacked vector (H), the entries 


Fi(yj,) = f f(x, y^) dx = 0 v27) 

XI 

for all i ® 1 , . . . , M. 

Having derived the matrix [T] which transforms the column vector 
{f} of data into the column vector (H), it is now a relatively easy 
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matter to derive an NM by NM matrix ISJ which gives a column vector of 
double integrals, that is. 


{A} » [S] {H} 


(28) 


where the NM by 1 column vector {A} is a stack of N, M by 1 column 
vectors {A } (m = 1 N) and 


f(X,y) dx dy, ..., 

yi *1 

f^M r*® f(x,y) dx dy)^ = (0, A 2 n„ Aj^, 
yi *^i 

T 


(29) 


In particular, if [J] is the M by M one-dimensional integrating matrix 
based on the grid G^, then 

{A^} = [J] (30) 


Equations (28) and (30) now show that [S] will be a banded matrix with 
N integrating matrices [J] along its diagonal, i-e., if [0] is an 
M by M zero matrix, then 


[S] 


[J] 

[0] 

N 

[0] 

[J] 

N 

\ 

. 

V 

X 




[ 0 ] [ 0 ] 


[01 


'[0] 

'tJ] 


(31) 


Combining equations (28) and (22) now shows that the required two-dimensional 
interating matrix [K], such that 


(A) = [K] {f} 


(32) 
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is the NM by NM matrix 


[K] = [S] [T] (33) 

As [S] and [T] involve only [Jl and the row vectors of [I], it is 
clear that [K] will be indopendent of data defining the function 
f(x,y) at the grid points. 

Sample matrices for [s], [t], and [k] associated with a 4 by 3 rectan- 
gular grid are given in appendix A. 

CONCLUDING REMARKS 


This report documents the derivation of two-dimensional integrating 
matrices for approximating integrals of a function of two variables on a 
rectangular region. The analysis shows that data defining the function 
at points on a rectangular grid which corresponds to the region must be used 
in a **stacked** column vector format. An appropriate two-dimensional integrating 
matrix is then obtained as the product of two auxiliary matrices which are 
formed from the one-dimensional integrating matrices corresponding to the two 
orthogonal directions on the grid* 
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APPENDIX A 


This appendix contains sample matrices [I], [J], 

[K] based on a 2-dimensional, rectangular grid G with 
(MM = 12, N = 4, M = 3). The x-grid G and y-grid G^ 

X V 

grid G are 


[S] , [T] , and 

12 points 
which define the 


G^ = {0., 12., 24., 36.} 


and 


= ^ 0 . , 2. , 4. } 

In these examples, the four by four integrating matrix [I] on the one- 
dimensional grid is based on vjuadratic interpolating pol>Tiomial s , while 

the three by three integrating matrix [J] on the one-Jimensional grid 
G^. is based on linear interpolating poUTioniial s . 
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Table A1 . 


The four by four integrating matrix [I] on the one -dimensional 

grid G based on quadratic interpolating polynomials. 

X 


. . 9 ^^ __CL^0 a.a. 

e.0 

10, a 16,0 0,0. 

- IB. a lSi. 0 ^6.0 5,0. 
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Table A3. The three by three integrating matrix [J] on the one-dimensional 
grid G based on linear interpolating polynomials. 


0.0 

170 
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• 4 dt 


Table A5 


The 12 by 12 integrating matrix [K] based on the two-dimensional 
rectangular grid G. 




0.0 

0.0 . 

0.0 

0,0 

9 , a 

0.0 

9,9 

0.0 

A.0 

0.9 


0 .^ 

0.0 

“«.0 

0.0 

' 0.0 

«.0 

0.0 

0.0 

0.0 

0.9 

0,0 


0«0 

. 0 * 0 .^ 


_ e * 0 _ 


9 . 0 - 

0,0 . 

0.9 

9.0 

0.9 

0,0 


0,0 

0.0 

0,0 

0.0 

0.0 

9 .«> 

0.0 

0.0 

0.0 

0,9 

0.0 



• 1.0 

~ 0,0 

5.0 

8,0 

- 1.9 

0,0 

9.9 

9.9 

0.0 

0.0 


6*0- 

-U0 



. 16»0 

• 3.0 

0.0 

5.0 

9.0 

• i.a 

0,0 


.0 

0.0 

0,0 

0,0 

0,0 

0.0 

0.0 

0.9 

9.9 

0.9 

0,0 


16«0 

• 2,0 

0,0 

10 .^ 

16.0 

• 3.0 

0.0 

0.0 

9.9 

0.0 

0,0 


lt.0 

- a .0 . 


.10 


• 4,0 

0,0 

19.9 

15.0 

•2.9 

0.0 


0.0 

0.0 

0.e 

0.0 

0 , J 

0.0 

0.0 

9.9 

0.0 

0.0 

0,0 

1<4«0 

15.0 

b.w 

” 5,0 

10,0 

15,0 

b .0 

5.0 

0.0 

0.9 

0,0 

0.0 


' 5«0 

i ^0 . 

S .0 

10 , 0 . 

30.0 

13.0 

10,0 

10.0 

is.e 

6.9 

5.0 
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